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We examine the isotope effect of superconductivity in systems with coexisting interactions 
of phonon and nonphonon mechanisms in addition to the direct Coulomb interaction. The 
interaction mediated by the spin fluctuations is discussed as an example of the nonphonon 
interaction. Extended formulas for the transition temperature T c and the isotope-effect coeffi- 
cient a are derived for cases (a) u np < u>o and (b) uj np > Wd, where u> np is an effective cutoff 
frequency of the nonphonon interaction that corresponds to the Debye frequency uj-q in the 
phonon interaction. In case (a), it is found that the nonphonon interaction does not change 
the condition for the inverse isotope effect, i.e., fi* > A p h/2, but it modifies the magnitude of 
a markedly. In particular, it is found that a giant isotope shift occurs when the phonon and 
nonphonon interactions cancel each other largely. For instance, strong critical spin fluctuations 
may give rise to the giant isotope effect. In case (b), it is found that the inverse isotope effect 
occurs only when the nonphonon interaction and the repulsive Coulomb interaction, in total 
effect, work as repulsive interactions against the superconductivity. We discuss the relevance of 
the present result to some organic superconductors, such as k-(ET)2Cu(NCS)2 and Sr2Ru04 
superconductors, in which inverse isotope effects have been observed, and briefly to high-T c 
cuprates, in which giant isotope effects have been observed. 
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The mechanisms of the pairing interactions of non- 
phonon origins have been examined for anisotropic su- 
perconductivity. For example, spin fluctuation exchange 
interactions have been examined for the superfluid 3 He 1 ^ 
and exotic superconductors, such as heavy- fermion su- 
perconductors, 2 - 1 organic superconductors, 3-5 ' high-T c 
cuprates, 6 ) and Sr 2 RuC>4 superconductors. 7 ) In order to 
clarify the origin of the pairing interactions, the isotope 
effect has been observed in many superconductors. It is 
widely known that the isotope-effect coefficient a = 0.5 
observed in nontransition-metal superconductors, such 
as Hg and Sn, is evidence of phonon-mediated pairing 
interactions in those systems. Here, the isotope-effect co- 
efficient a has been defined by 

dlnT c 

where M is the relevant atomic mass. The value 
a = 0.5 is directly derived from the BCS formula 
T c = 1.13w D e- 1 / A with w D oc Af" 1 / 2 and A oc M°. 

However, deviation of a from 0.5 is not necessarily 
evidence of the presence of the pairing interaction of a 
nonphonon origin. For instance, in transition-metal su- 
perconductors, such as Ru and Os, the value of a largely 
deviates from 0.5 due to the strong repulsive Coulomb 
interaction Vq- In the presence of it, A is replaced with 
A = A — /j* , where 



Hence, we obtain 



l + /iln(M//w D ) : 



(2) 



/i = VcN(0) and W denotes an effective cutoff energy of 
the Coulomb interactions of the order of the band width. 



1 



(3) 



In particular, the strong Coulomb repulsion such as 
A > [i* > A/2 gives rise to the inverse isotope effect 
(a < 0). The anharmonicity of the lattice vibrations 
may also affect the isotope effect in compounds, such 
as PdH(D) (a « — 0.25). 8 ) The inverse isotope effect has 
been observed also in some organic superconductors 9-15 ^ 
and Sr2Ru04, 16 ) but the mechanisms is unknown. 

In most cases, a nonzero isotope shift (a ^ 0) suggests 
the presence of the phonon contribution to the pairing 
interactions. One might expect that isotope shift can be 
attributed to the change in normal state properties, such 
as the density of states. However, usually, such changes 
should be of the order of (AM/2M)(w D /e F ), which is 
usually negligible except for very narrow band systems, 
where AM is the shift of the atomic mass. Therefore, 
the isotope effect is mainly attributed to a change in the 
pairing interactions. 

Here, it could not be excluded a priori that a ^ oc- 
curs because the phonon interaction Vph(p,p') reduces 
the anisotropic pairing interaction, but it seems unlikely. 
The effective interaction for the anisotropic gap func- 
tion proportional to 7(p) is essentially proportional to 
the average of V p h(p, p') 7*(p)7(p') on the Fermi sur- 
face, where 7(p) is a function for expressing the momen- 
tum dependence. Since Vph(p,p') tends to take larger 
negative values for smaller momentum transfers, i.e., for 
p ps p', where 7*(p)7(p') > 0, the above average must 
be negative. For example, in order to obtain a repulsive 
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effective interaction for 7(p) oc p x , Vph(p,p') must have 
a negative peak at |p — p'| ~ 2fcp rather than |p — p'| ^ 0, 
but it is unlikely. In fact, some explicit calculations have 
shown that the phonon interaction includes attractive 
components for anisotropic pairing. 17 ~ 19 ) Therefore, the 
magnitude of the isotope shift of T c reflects, more or less, 
the degree of the positive contribution to the Cooper pair 
formation from the phonons. 

We should note that spin triplet pairing does not nec- 
essarily exclude the phonon mechanisms of the pair- 
ing interaction. The triplet pairing interaction is usually 
weaker than the singlet pairing interaction in the phonon 
mechanism, but it could become dominant, for example, 
in the presence of a strong short-range Coulomb interac- 
tion, 17 ~ 19 ) which strongly suppresses the isotropic pairing 
more than the anisotropic pairing. In this context, the 
large deviation of a from 0.5 is a natural consequence of 
the strong Coulomb repulsion. The layer structure also 
favors the anisotropic pairing interactions mediated by 
phonons. 19 ) 

In this paper, we examine the isotope effect in systems 
in which interactions of phonon and nonphonon mecha- 
nisms coexist. In particular, we examine the condition 
for the inverse isotope effect. 

We assume that the nonphonon interaction also has 
an effective cutoff frequency uj np that corresponds to 
the Debye frequency u>d in the phonon interaction. This 
assumption is explained as follows for the interactions 
mediated by the spin fluctuations. The spin fluctua- 
tions have a characteristic frequency u s f , which becomes 
smaller as one approaches the magnetic instability point 
for the critical slowing down. The structure of the suscep- 
tibility x(q, uo) reflects this characteristic frequency. The 
susceptibility x(q, has a sharp peak around a momen- 
tum q with a width corresponding to w s f, in addition 
to a broad peak over the whole momentum space, which 
reflects the electron characteristic energy ep. Since the 
effective interaction mediated by the spin fluctuations 
is essentially proportional to x(qjw), it consists of two 
parts with different characteristic energies u) s { and ep • In 
the gap equation, these characteristic energies play roles 
of the effective cutoff energies of the two parts of the ef- 
fective interaction. The part of the effective interaction 
with ep can be included in the repulsive Coulomb inter- 
action, which also has a large effective cutoff energy of 
the same order. Therefore, we may regard that the inter- 
action mediated by the spin fluctuations has an effective 
cutoff energy of the order of ui s f, namely, uj np ~ u s f. 
In fact, it has been shown in a theoretical calculation 5 -* 
that the contribution to the pairing interaction mainly 
comes from this sharp peak of x- 111 theoretical phase 
diagrams, 5 ' 20 -' the superconductivity by such a pairing 
interaction occurs only near the magnetic phase bound- 
ary. An explicit calculation, which illustrates the above 
argument, has been given in ref. 5. 

We start with the interactions between electrons, 



coordinate on the Fermi surface and the electron energy 
measured from the Fermi surface with indices k and a 
to express the kinds of interactions and the pairing sym- 
metries. The symmetry functions J a (p\\) are normalized 
with respect to the average on the Fermi surface. 21 ) We 
consider a gap function proportional to 7 Q (p||), retaining 
only terms with a in eq. (4). Introducing cutoff energies 
LJk, we put 



(5) 



where io\ < lo-i < u>3 < ■ ■ ■ < uj n . We should note that 
V k a n for the repulsive Coulomb interaction can be nonzero 
even for anisotropic pairing, unless we assume only on- 
site repulsion. For eq. (5), the gap function is written 

as 



A a (p) = A k» 6 ( UJ k - |£p|)7a(P||)- 



(6) 



fe=i 



The linearized gap equation (T c equation) is written as 

k—l n 
Afe Q = Afc|^^ &-k'uh' + ^2 A k'alk , (7) 
k' = l k'=k 

where we have defined A& = — ly^7V(0) and l k — 
ln(2e 7 o) fe /7rT c ) with the Eulcr constant 7 = 0.57721 
We write the transition temperature in the form 



2eT 

T c = u)ie 

TT 



1/Ai 



(8) 



defining the effective coupling constant Ai. We define 
ako = —dlnujk/dlnM for convenience. 

First, we consider the case with n = 2. The T c equation 
is written as 



Ar 1 ~h 

-h 



-h 

A2 - h 



Ala 
A 2 « 



(9) 



Therefore, we obtain Ai = Ai + A2 and A2 = A2/[l — 
The isotope-effect coefficient is obtained 



A 2 ln(u> 2 /cji) 
as 



a 



aw 



AS 



Ai + A* 



«20 



Ai + Ai 



(10) 



v( P ,p) = }^}^v k a (^c P 'h a (p\\h a (p\\ ), (4) For case (a); we put Wl 



If we take and V20 as a phonon-mediated pairing in- 
teraction and the repulsive Coulomb interaction, respec- 
tively, then oj\ — ojv oc M -1 / 2 and w 2 = W oc M°, and 
eq. (10) is reduced to the standard formula (3). 

Equation (10) can be used for the systems in which two 
different interactions coexist and the repulsive Coulomb 
interaction is negligible. Here, we consider a system in 
which one of the two interactions is of a nonphonon 
origin. We examine two cases: (a) uj np < wd and (b) 
wd < w np , where u> np oc M° and wd oc M- 1 ' 2 . If the 
nonphonon interaction is one mediated by the spin fluc- 
tuations, uj np ~ lu s {, and case (a) occurs in proximity to 
the magnetic instability. We define A p h = —V°^N(0) and 
A np = — V" p N(0), where V p \ and V" p denote the coupling 
constants of the phonon and nonphonon interactions, re- 
spectively. 



a k=l 



-'npj 



UJ 2 



wd, Ai — A 



np i 



where p\\ and £p denote the two-dimensional momentum 
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Fig. 1. Isotope-effect coefficients as functions of A p h for case (a) 
LU np < wrj- Solid, dashed, and dotted curves show the resuls for 
A np = 0.2, 0, and —0.2, respectively. Thick and thin curves show 
the values of a and Ai, respectively. 



Fig. 2. Phase diagram for case (b) L0 np > u>d, with oj np /wrj = 2 
and fi* = 0.1. NSC denotes the region where the superconduc- 
tivity does not occur. 



and AS = A* h . Thus, we obtain 



1 



a 



A P h 



(11) 



2 A np + A* h ' 

which is positive definite. On the other hand, for case 
(b), we put loi = w D , w 2 = w np , Ai = A ph , and A 2 = A np . 
Thus, we obtain 

a; 



a= 2 



l-( 



'np 



Aph + A* F 



(12) 



where A* p = A np /[1 - A np ln(w np /a) D )]. Therefore, a 
could become negative only in case (b) (uj np > cjd) with 
A p hA* p < 0. 

Next, we examine the effect of the direct Coulomb in- 
teraction. The T c equation with n — 3 can be solved in a 
similar manner to that in the above. It is obtained that 
Ai = Ai + A 2 , AS = (A 2 + - (A 2 + A3) hi], and 

A3 = As/[1 — A3/3 2 ], where Ikw = l n ( a; fc/<- t 'fc)- The transi- 
tion temperature is given by eq. (8). Hence, we obtain 



a = a w Ci + a 2 oC 2 + a 30 C 3 , 



(13) 



where 



Ci = 1 - ( 



AX 



C 2 = 



c 3 = ( 



y 2 



\l + \* 2 

2 



Ai + A* 

A 2 



A 2 + A* 



(14) 



) 2 (- 



A 3 



Ai + AS A 2 + A 



We examine cases (a) and (b) mentioned above, adding 
the direct Coulomb interaction Vg" with lj 3 — W and 
A 3 = -fi. 



In case (a), we obtain Ai = A np + A* h and 



2^A 



A 



ph 



) 2 i-( 



A p h - A* 



(15) 



where 
with \i 



A ph = (\h - M*)/[l - (Aph - H*) ln(w D /w np )] 
denned by eq. (2). In this case, the sign of a 
is determined only by the ratio A p h/yu*- If we assume 
fi* > 0, there is no region of A p h < where a < 



occurs. In Fig. 1 and eq. (15), it is also found that 
a rapidly varies with A p h and A np . The prefactor 
[A p h/( A np + Aph)] 2 m C( l- (15) gives rise to a large value 
of a when A* h w — A np . 

On the other hand, in case (b), we obtain Ai = A p h + 
A* p and 



a= 2 



l-( 



A* 



Aph + A* p ' 



(16) 



where A* p = (A np - - (A np - //*) ln(w np /w D )] 

and in this case [i* = + fi\n(W/uj n p)]. It is found 
that A p hA* p < is a necessary condition for a < 0. 
Since it is reasonable to assume A p h > even for 
anisotropic pairing as argued above, an inverse isotope 
effect suggests that the nonphonon interaction and 
the direct Coulomb interaction, in total effect, work 
as repulsive interactions against the superconductivity 
(A* p < 0) when uj np > uj^,. Figure 2 shows the phase 
diagram obtained using eq. (16). 

Now, we describe the general solution of eq. (7) with 
any positive integer n. Equation (7) is rewritten as = 
AfeFfe, where we have defined = X)fc'=i h'^k' + h^-k 
with Afc = X]fe'=fe^fe'- On the other hand, it is easily 
proved from eq. (7) with k + 1 <~ n that an equation 
of the form Afe + i = A£ +1 i<fc holds for k > 1 with a 
recurrence formula 



A k — 



A fc + At 



fc+i 



1 — (Afe + *l +1 )lk,k- 



(17) 



for k ^ 2. Hence, it holds that Ai = Ai + A 2 = (Ai + 
A^ZiAi. Therefore, T c is given by eq. (8) with Ai = 
Ai + Aj. 22 ) Hence, from eq. (1), the general formula for 
a is obtained as 



a 



(18) 



k=i 



with Ck — Afe — Afc + i for k < n — 1 and C n 

jk — 1 f\* 1 1 \ _ 1 \ * M 2 



A„, where 
1. 



Afc = ^-iZiiK+ili^i + Af+i)] 2 for k > 2 and Ai 
For example, if the interactions V£ with k — 1, ■ • -n — 1 
are phonon interactions with otko = 1/2, it is found that 
a w 1/2 for large n, since A„<1. 
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Let us discuss the inverse isotope effect observed in 
clean samples of Sr 2 Ru0 4 with T c > 0.94T c0 , 16) where 
T c0 is the transition temperature in the clean limit. From 
the present result for case (b), which is probably the case 
for Sr 2 RuC>4, 23 ) the inverse isotope effect suggests that 
the superconductivity in this compound is induced by 
the phonon-mediated pairing interaction. 24 ' 

This result is consistent with the recent report that 
the antiferromagnetic fluctuation does not seem to play 
a constructive role in the triplet superconductivity in 
Sr 2 RuC>4. 25 ) The impurity concentration dependence of 
the isotope-effect coefficient 16 ' can also be reproduced 
within a theory based on a phonon mechanism. 21 - 1 Many 
novel mechanisms have been proposed for the supercon- 
ductivity in this compound. 26 ' However, there is no con- 
clusive evidence for the nonphonon mechanism at the 
present. The present theory does not consider the effect 
of anharmonicity of lattice vibrations. However, if such 
effect is essential for the inverse isotope effect, we also 
obtain the same conclusion that the phonon interaction 
is dominant. 

The situation is rather different in the organic 
superconductors. The inverse isotope effect was 
also observed in the organic superconductors, such 
as K -(ET) 2 Cu(NCS) 2 , 10 ' 11 ' 13 ) /3-(ET) 2 I 3 , n ' 12 ' and 
kl- (ET) 2 Ag(CF 3 ) 4 ( 1-bromo-f ,2-dichloroethane) , 9 ' 
while a normal isotope effect was observed in k- 
(ET) 2 Cu[N(CN) 2 ]Br n < 14 < 15 ) In contrast to Sr 2 Ru0 4 , 
case (a) (u s f < ujd < W) may be realized in the organic 
superconductors even if £ mag = 2a ~ 3a, because they 
have narrower band widths of 0.2eV ~ 0.5cV. The 
coefficient a varies rapidly, as shown in Fig. 1, when A p h 
changes, whichever A np > or < 0. Thus, possibilities of 
the nonphonon pairing interactions arc not excluded in 
the organic superconductors with a < 0. 

In high-T c cuprates La 2 _ rE Sr :E Cu04, it was observed 
by Crawford et aZ. 27 ' that a increases near the antifer- 
romagnetic boundary, where T c decreases. In cuprates, 
the presence of the strong spin fluctuations has often 
been pointed out especially for underdoped samples. For 
example, we have discussed a possible relevance of the 
strong spin fluctuation to the suppression of T c in an 
underdoped region. 28 ' Therefore, the giant isotope shift 
in the present result for case (a) may explain the experi- 
mental results. A more detailed discussion on this subject 
would be presented in a separate paper. 

In conclusion, we have proposed a model of supercon- 
ductors with coexisting interactions of phonon and non- 
phonon origins. In the presence of the Coulomb interac- 
tion, the extended formulas for T c and a, and a phase 
diagram are obtained. It is found that the inverse iso- 
tope effect occurs due to the cancellation of the attrac- 
tive and repulsive interactions, where one of them needs 
to be the phonon interaction, while the other the non- 
phonon interaction with a larger effective cutoff energy 
than the Debye frequency. When the phonon interaction 
is attractive, any additional nonphonon pairing interac- 
tion reduces the value of a, but does not make a negative. 
It is also found that a strong critical spin fluctuation with 
oj s f < ujd could give rise to a giant isotope shift. 
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